Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



XMM- Newton VIEW OF THE ULTRA-LUMINOUS X-RAY SOURCES IN M 51 
GuLAB C. Dewangan\ Richard E. Griffiths\ Manojendu Choudhury^, Takamitsu Miyaji^, & Nicholas J. 

SCHURCH^ 
Draft version February 5, 2008 

ABSTRACT 

We present results based on XMM-Newton observation of the nearby spiral galaxy M 51 (NGC 5194 
and NGC 5195). We confirm the presence of the seven known ULXs with luminosities exceeding the 
Eddington luminosity for a 1OM0 black hole, a low luminosity active galactic nucleus with 2 — 10 keV 
luminosity of 1.6 x 10'^^ erg s~^, and soft thermal extended emission from NGC 5194 detected with 
Chandra. In addition, we also detected a new ULX with luminosity ^ 10^^ erg s~^. We have studied 
the spectral and temporal properties of the LLAGN and 8 ULXs in NGC 5194, and an ULX in 
NGC 5195. Two ULXs in NGC 5194 show evidence for short-term variability, and all but two ULXs 
vary on long time scales (over a baseline of ~ 2.5 years), providing strong evidence that these are 
accreting sources. One ULX in NGC 5194, source 69, shows possible periodic behavior in its X-ray 
flux. We derive a period of 5925 ± 200 s at a confidence level of 95%, based on three cycles. This 
period is lower than the period of 7620 ± 500 s derived from a Chandra observation in 2000. The 
higher effective area of XMM-Newton enables us to identify multiple components in the spectra of 
ULXs. Most ULXs require at least two components - a power law and a soft X-ray excess component 
which is modeled by an optically thin plasma or multicolor disk blackbody (MCD). However, the 
soft excess emission, inferred from all ULXs except source 69, are unlikely to be physically associated 
with the ULXs as their strengths are comparable to that of the surrounding diffuse emission. The 
soft excess emission of source 69 is well described either by a two temperature mekal plasma or a 
single temperature mekal plasma {kT ~ 690 eV) and an MCD {kT ~ 170 eV). The MCD component 
suggests a cooler accretion disks compared to that in Galactic X-ray binaries and consistent with 
that expected for intermediate mass black holes (IMBHs). An iron Ka line {EW ^ 700 eV) or K 
absorption edge at 7.1 keV is present in the EPIC PN spectrum of source 26. The spectrum of the 
ULX in NGC 5195, source 12, is consistent with a simple power law. The LLAGN in NGC 5194 shows 
an extremely flat hard X-ray power-law (F ~ 0.7), a narrow iron Ka line at 6.4 keV (EW ^ 3 keV), 
and strong soft X-ray excess emission. The full band spectrum is well described by a two component 
mekal plasma and reflection from cold material such as putative torus. 

Subject headings: accretion, accretion disks — galaxies: active — galaxies: individual (M 51) — 
X-rays: binaries — X-rays: galaxies 



1. INTRODUCTION 

X-ray emission from spiral galaxies consists of sev- 
eral components including an active nucleus (if present), 
extra-nuclear point X-ray sources, hot gas, supernova etc. 
Many spiral galaxies contain a mildly or non-active nu- 
cleus and, in this case, the luminosity in the X-ray band 
is dominated by the emission from a few luminous off- 
nuclear point sources (Roberts & Warwick 2000). The 
brightest of these X-ray sources, with luminosities ex- 
ceeding the Eddington luminosity for a 1OM0 black hole, 
are defined as the ultra- luminous X-ray sources (ULXs). 
These objects, which may sometimes individually out- 
shine the rest of the galaxy, were first discovered in Ein- 
stein observations (Fabbiano & Trincheri 1987; Fabbiano 
1989). Later they were detected in large numbers by 
ROSAT (Colbert & Mushotzky 1999; Roberts & War- 
wick 2000; Colbert & Ptak 2002). The first detailed 
X-ray spectral characteristics of ULXs were revealed by 
ASCA and demonstrated that most ULXs exhibit mul- 
ticolor disk blackbody spectra with a 'high temperature' 
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similar to that observed from Galactic X-ray binaries in 
their high state. Some ULXs were also found to vary 
and make a transition, at lower luminosities, to power- 
law spectra, similar to the spectral transitions seen in 
many Galactic XRBs. However, the observed X-ray 
flux of many ULXs implies (assuming isotropic emission 
and Eddington-luminosity limited systems) black holes 
of ~ 10^ — 1O^M0. These masses are much larger than 
the measured masses of known Galactic X-ray binaries 
or black holes (GXRBs or GBHs) and much lower than 
the measured masses of active galactic nuclei, thus mak- 
ing them candidates for intermediate-mass black holes 
(IMBHs). 

Several models have been proposed to explain the high 
luminosities of ULXs. These include {i) accreting IMBHs 
with sub-Eddington rates (e.g., Colbert & Mushotzky 
1999), (ii) XRBs with anisotropic emission (King et al. 
2001), {Hi) beamed XRBs with relativistic jets directly 
pointing towards us i. e., scaled down versions of blazars 
(Mirabel & Rodriguez 1999), and (iv) XRBs with super- 
Eddington accretion rates (Begelman 2002). A few ULXs 
have been found to be young supernova (e.g., SN 1979, 
Immler, Pietsch & Aschenbach 1998) with X-ray emis- 
sion powered by the interactions of shock waves and the 
surrounding dense inter stellar medium. Some of the 
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accretion-powered ULXs are also observed to be associ- 
ated with older supernova remnant-like structures (e.g., 
IC 342-Xl, Roberts et al. 2003). The recent observa- 
tion of a break at a frequency of ~ 28 mHz in the power 
density spectrum of the ULX NGC 4559 X-7 suggests a 
mass of a few thousands solar masses, consistent with the 
measured thermal temperature of fcT — 120 eV. This ob- 
servation clearly supports the IMBH scenario for ULXs. 
It has become increasingly clear that there are two sepa- 
rate populations of ULXs. A large number of ULXs are 
located in regions of current star formation in star-burst 
and spiral galaxies (The Antennae, Zezas et al. 2002; 
NGC 3256, Lira et al. 2002; NGC 4485/90, Roberts et 
al. 2002; Arp 299, Zezas, Ward, & Murray 2003; M 51, 
Terashima & Wilson 2004). These ULXs are likely to be 
associated with young stellar populations. Other ULXs 
have been observed in early- type elliptical/SO galaxies 
and are likely to be associated with older stellar pop- 
ulations. A recent study of the ULX population of 82 
galaxies by Swartz et al. (2004) suggests that ^ 25% of 
these sources may be background objects, including 14% 
in spiral galaxies and 44% in elliptical galaxies. 

Chandra has revolutionized the detection of ULXs pri- 
marily due to its superior angular resolution (see e. g., 
Swartz et al. 2004). XMM-Newton, on the other hand, 
has been probing the detailed characteristics of the ULX 
populations in nearby galaxies due to its large effec- 
tive area. XMM-Newton is well suited to study detailed 
spectral and temporal characteristics of ULXs in nearby 
galaxies (e.g., Strohmayer & Mushotzky 2003; Cropper 
et al. 2003; Miller et al. 2004; Dewangan et al. 2004). 

M 51, also known as the Whirlpool galaxy, is a face- 
on spiral galaxy and is located at a distance of 8.4 Mpc 
(Feldmeier et al. 1997). Optical studies of emission lines 
classified it as a LINER or a Seyfert 2 galaxy (Stauffer 
1982). Furthermore, Ho et al. (1997) suggest the pres- 
ence of a broad Ha line. Kohno et al. (1996) found a 
nuclear molecular disk, and constrained the dynamical 
mass within 70 pc of the nucleus to be 4 — 7 x 10^ Mq. 
This indicates that M 51 hosts as massive a black hole 
as many AGNs. 

X-ray observations also strongly suggest the presence 
of an AGN in M 51. Einstein (Palumbo et al. 1985) and 
ROSAT (Marston et al. 1995; Ehle et al. 1995) observa- 
tions constrained the soft X-ray luminosity of the M 51 
nucleus to be Lx < 5 x lO^^erg s~^. Ginga scanning 
observations in 1988 detected bright hard X-ray emis- 
sion with 2 — 20 keV luminosity of 1.2 x 10^^ erg s~^, 
a photon index of 1.43 ± 0.08 and an intrinsic absorp- 
tion of < 7 X lO^^cm-^ (Makishima et al. 1990). Such 
hard nuclear X-ray emission is often considered to be 
evidence for a low luminosity AGN (LLAGN). ASCA 
observed M 51 in the hard (2 — 10 keV) X-ray band 
in 1993, and did not detect such a bright hard com- 
ponent. Instead, a faint hard X-ray continuum with a 
neutral Fe Ka line was detected, whose flux was an or- 
der of magnitude lower than that measured by Ginga 
(Terashima et al. 1998). The reason for this large dif- 
ference in luminosity has recently been clarified by Bep- 
poSAX observations. Fukazawa et al. (2001) find that 
the nucleus is photoelectrically absorbed below 10 keV 
but is observed directly above 20 keV, implying an ab- 
sorbing column of Nh — 5.6 x lO^'* cm~^. The 2-10 keV 
luminosity measured by BeppoSAX was similar to that 



of ASCA. Fukazawa et al. (2001) attributed the higher 
Ginga luminosity to variability of the absorbing column. 

High angular resolution ROSATHYd observations also 
revealed eight X-ray point sources and diffuse soft X-ray 
emission in M 51. Chandra observed M 51 during June 
2000 and June 2001 with ACIS-S. These data have been 
studied in detail by Terashima & Wilson (2001, 2004). 
The X-ray image revealed the nucleus, 113 X-ray sources 
and extended emission which resembled the morphol- 
ogy of both the radio and optical emission line images 
(Terashima & Wilson 2001). The X-ray image of the 
nucleus is well represented by a model consisting of soft 
thermal plasma {kT ~ 0.5 keV), a very hard continuum 
and an Fe Ka emission line at 6.45 keV with an equiv- 
alent width of greater than 2 keV. The X-ray spectra 
of the extra-nuclear clouds are well fitted by a thermal 
plasma model with kT ~ 0.5 keV. The spectral shape 
and morphology strongly suggest that the clouds are 
shock-heated by a bipolar outflow from the nucleus. Out 
of 113 extra- nuclear sources, 9 sources have luminosities 
exceeding 10'^^ erg s~^ in the 0.5-8 keV band. The num- 
ber of extra- nuclear sources in M 51 is much higher than 
in most other normal spiral and elliptical galaxies and is 
similar to galaxies experiencing star-burst activity. The 
X-ray spectra of most of the detected sources are consis- 
tent with a power-law spectral form with a photon index 
between 1 and 2, while one source has an extremely hard 
spectrum and two sources have particularly soft spectra. 
The X-ray spectra of three of the ULXs are consistent 
with both a power-law and multi-color blackbody model. 
One ULX showed a remarkable spectrum with prominent 
emission lines. 

In this paper we present an XMM-Newton EPIC ob- 
servation of M 51. Utilizing the large collecting area of 
XMM-Newton, we have studied spectral and temporal 
properties of nine ULXs and the LLAGN in M 51. The 
paper is organized as follows. In Sect. 2, we outline the 
observation and the data selection. In Sect. 3, 4, and 
5, we present the spatial, temporal, spectral analysis. 
Sect. 6 describes properties of the individual ULXs. We 
discuss the results in Sect. 7, followed by conclusions in 
Sect. 8. 

2. OBSERVATIONS AND DATA SELECTION 

The galaxy M 51 was observed by XMM-Newton ob- 
servatory ( Janscn ct al. 2001) on 2003 Jajnuary 15 for a 
duration o f ~ 21 ks. The EPIC PN ifStru'der et al. 200Hl 
and MOS l|Turner et al. 200 1) cameras were operated in 
full-frame mode using the thin filter. The optical monitor 
(OM) instrument and the reflection grating spectrome- 
ters operated simultaneously with the EPIC cameras. 

The raw PN and MOS events were processed and fll- 
tered using the most recent version of the calibration 
database and analysis software (SAS v6 . 0) available in 
March 2004. Examination of the background light curves 
extracted from source-free regions showed that our ob- 
servations was not strongly affected by particle induced 
flares. Events in known bad pixels were discarded. For 
the spectral analysis, we used events with patterns 0-4 
(single and double) and flag zero for the PN, and pat- 
terns 0-12 (similar to ASCA event grades 0-4) and flag 
zero for the MOS. However, for the temporal analysis, 
we used all events with pattern 0-12 for both the PN 
and MOS cameras. 
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3. SPATIAL ANALYSIS 

We extracted PN images in the full (0.2 - 10 keV), 
0.2 - 0.7 keV, 0.7 - 2 keV, and 2-10 keV bands. First 
we adaptively smoothed the full band image using the 
CIAO task csmooth and created a smoothing kernel. The 
images in the 0.2 - 0.7 keV, 0.7 - 2 keV and 2-10 keV 
bands were then adaptively smoothed using the kernel 
created above, and combined to form a three-color im- 
age shown in Figure ^ (Left). The color image shows 
several discrete sources, diffuse soft emission, and the 
bright nucleus of NGC 5194. Interestingly, the bright 
nucleus of NGC 5194 actually consists of two hard X- 
ray sources ~ 20" apart (see below). Also present is the 
companion galaxy NGC 5195 showing diffuse emission, 
multiple sources in the central region, and a bright point 
source to the left of the nucleus. To avoid any possible 
confusion on multiple source identification numbers aris- 
ing from multiple X-ray missions, we follow the Chandra 
source identification numbers as set by Terashima & Wil- 
son (2004). Figure^{right) shows the EPIC PN image in 
the 2 — 10 keV band and the positions of the ULXs stud- 
ied in this paper have been marked. The ULX source 9 in 
NGC 5194 is not seen in the hard band image image due 
to its location in the chip gap and its extreme softness. 
The source counts and positions listed in Table 1 have 
been determined by a source detection procedure by uti- 
lizing a series of SAS tasks as described by Miyaji et al. 
(2003). In short, an input image is searched for point- 
like sources by a simple sliding-cell detection method 
with local background using eboxdetect. This is followed 
by a "background map" creation by a spline fit to the 
source-excluded region using esplinemap. A sliding cell 
detection is repeated to find excess over this background 
map (eboxdetect). The final source positions and counts 
are determined by a multi-source maximum-likelihood fit 
with the XMM-Newton PSF, using the sources detected 
in the previous step as the starting points (emldetect). 
In order to be free from the effects of soft diffuse emis- 
sion, which could prevent the accurate determination of 
the background map, the hard-band (2-10 keV) image is 
used to obtain the source counts, except source 9, which 
has not been detected in the hard band. In Figure El 
we show the UV image of M 51 taken with the XMM- 
Newton OM. The image shows beautiful spiral arms, a 
prominent nucleus, and bright star clusters surrounding 
the nucleus with either a ring-like morphology or multi- 
ple inner spiral arms. One of the spiral arm of NGC 5194 
appears to be tidally connected to the companion galaxy 
NGC 5195. This galaxy shows a bright nucleus with dark 
dust lanes near the nucleus. 

In Fig. 121 we have overlaid the X-ray contours onto 
the UV image. The contours represent the hard (1.5 — 
12 keV) X-ray sources only. Fig. |21 highlights the 
two bright hard X-ray sources in the central region of 
NGC 5194. The two sources are separated by ~ 20". 
The fainter of the two X-ray sources coincides with the 
UV nucleus. The brighter source is an ULX, which is 
clearly resolved in Chandra images (Terashima & Wil- 
son 2004). 

In this paper we study the active nucleus and the 
point sources with luminosities > lO'^^ erg s~^ in the 
0.3 — 12 keV band. There are eight ULXs in addition 
to the LLAGN in NGC 5194. There is also a well iso- 



lated ULX in the companion galaxy NGC 5195. How- 
ever, the central region of NGC 5195 consists of multiple 
bright point sources and extended soft emission as seen in 
the Chandra ACIS-S image (Terashima & Wilson 2004). 
Due to lower angular resolution of XMM-Newton, it is 
not possible to extract counts for individual sources in 
the central region of NGC 5195. Therefore we exclude 
the bright sources in the central region of NGC 5195 from 
our analysis. The Chandra source NGC 5194#63 had an 
X-ray luminosity < 10"^^ erg s~^ during the two previous 
Chandra observations. In this observation this source has 
a luminosity of ~ 2 x 10^^ erg s~^ in the 0.3 — 10 keV 
band and we now include this source in our sample. 

4. TEMPORAL ANALYSIS 
4.1. Short-term variability 

We extracted 0.2 — 12 keV X-ray light curves for the 
LLAGN and ULXs in M 51 from the EPIC PN data. The 
extraction regions had radii in the range of 15 — 35" and 
enclosed > 60% of the X-ray emission from each source. 
The smallest circular regions were chosen to avoid con- 
tamination from nearby sources. We used time bins of 
520 s and required that the bins be at least 50% exposed 
in order to minimize the errors due to the low number of 
counts from the weak sources. We also extracted individ- 
ual background light curves for each of the sources using 
circular regions lying approximately at the same readout 
distance as that for the source. The source light curves 
were then corrected for the background contribution us- 
ing the background light curves and appropriate scaling 
factors to take into account the different extraction areas. 

We fit a constant count rate model to each of the light 
curves to investigate the possibility of rapid variability. 
We rebin the light curves, if required, to ensure a mini- 
mum of 20 counts per bin to allow the use of fitting 
statistics. The results for each light curve are listed in 
Table 2. We find that only two of the sources show sig- 
nificant variability: source 26 (at the 98.6% confidence 
level) and source 69 (99.95% confidence level). The light 
curves of these two sources are plotted in Figure |31 

4.2. Periodicity 

Source 69 shows the most significant rapid variability. 
The light curve shows three tentative peaks at regular 
intervals. The first Chandra observation of this ULX 
in June 2000 showed two peaks separated by ~ 7000 s. 
These features were interpreted by Liu et al. (2002) as 
a signal with a period of ^ 2.1 hr. We have searched 
for a similar periodic signal in the XMM-Newton light 
curves. For this purpose, we combined the PN and MOS 
light curves, binned to 520 s. The initial search for pe- 
riodicity in the source, as well as the folding of the light 
curve with respect to the final period obtained was per- 
formed with the XRONOS 5.19 package, which is a part 
of FTOOLS 5.2. The periodogram was obtained using 
the algorithm of Scargle (1982), following the recipe of 
normalizing the periodogram power using the total vari- 
ance of the complete light curve as prescribed by Horne 
& Baliunas (1986). This algorithm is very robust for 
detecting periodic signals in light curves with few data 
points as well as non-uniformly distributed observational 
data, while providing a 'false alarm probability' of the 
periodicity arising in the light curve from random noise. 



4 



This probability is obtained using the signal to noise ra- 
tio, where the variance of the noise is obtained under 
the basic assumption that it is a normally distributed 
random variable. The periodogram is shown in Figure^ 
{left). The period obtained was 5925 ±200 s, with a 'false 
alarm probability' of 0.045, suggesting that the observed 
modulation is not generated by the statistical fluctuation 
of the noise, to a confidence level of 95.5%. The main 
source of the noise in the power in the Fourier frequency 
domain are the inherent statistical fluctuations at scales 
< 2 ks (Figure 01 . We have tested for the periodicity 
on top of a constant count rate. However, in reality, the 
source probably has a red-noise power spectrum, simi- 
lar to XRBs. Therefore, the ~ 2a or 95.5% detection of 
the period is an upper limit, and the periodicity may or 
may not be a true signal. Future long baseline X-ray ob- 
servations will be crucial for the robust detection of the 
periodicity from this ULX. In Figxn:e^{right) we provide 
the hght curve folded with respect to the period of 5925 
s, with a bin size of 520 s. Evidently, the modulation 
fraction is in excess of 45%, similar to the estimate of 
Liu et al. (2002). 

4.3. Long-term Variability 

We have compared the 0.5 — 8 keV fluxes derived from 
the XMM-Newton observation (see below) and earlier 
two observations with Chandra (Terashima & Wilson 
2004) in June 2000 and June 2001. Figure El shows 
the comparison of the fluxes. Cross-calibration results 
suggest that the fluxes measured with Chandra and 
XMM-Newton agree at a level better than 10% (Snow- 
den 2002). Thus, all the ULXs, except NGC 5194#41 
and NGC 5195#12, show long-term significant variabil- 
ity over a baseline of 2.5 yr. 

5. SPECTRAL ANALYSIS 

For each ULX and the LLAGN, we extracted source 
and background spectra from the final filtered PN event 
list using the same extraction regions that we used for 
the PN light curves. We also extracted the source spec- 
tra from the MOS data using source regions similar to 
that used for the PN. The background MOS spectra for 
the individual ULXs were extracted using circles in the 
nearby source-free regions. Appropriate PN and MOS 
responses and effective area files were created for each 
source using the SAS tasks rmfgen and arfgen. A group- 
ing of 20 counts per spectral channel resulted in signif- 
icant degradation in the spectral resolution of the ULX 
spectra above 2kev. Therefore, all the source spectra 
were grouped to a minimum of 10 counts per spectral 
channel so that discrete spectral features are not missed. 
The spectra were analyzed with the XSPEC 11.3 spec- 
tral analysis package (Arnaud 1996), and using the C- 
statistics instead of x^-statistic which is not useful in 
the case of low ( ^ 20) counts per bin. Unlike the x^- 
statistic, the C-statistic does not provide a goodness-of- 
fit (GOF) criterion. Therefore we simulate 5000 spectra 
based on the best-fit model derived using the C-statistic, 
and calculate the percentage of these spectra with the 
C-statistic less than that for the data. This number pro- 
vides the GOF, and if this number is ^ 50%, then the 
model is a good fit to the observed spectrum (similar to 
the reduced 1 criterion), though obviously we can- 
not statistically exclude fits for which GOF< 95%. The 



confidence intervals are calculated in the same way as in 
the case of x^-statistic. All the errors quoted below were 
calculated at 90% confidence level for one interesting pa- 
rameter i.e., AC = 2.7. 

At first the MOS and PN spectra of a few bright 
sources were fitted separately to check for the possi- 
ble uncertainties due to cross-calibration problems. We 
found generally good agreement between MOS and PN 
cameras in the 0.3 — 10 keV bands. Therefore we present 
the spectral results obtained by fitting the same model 
jointly to the PN and MOS data while leaving the rela- 
tive normalizations to vary. 

5.1. The low luminosity AGN in NCC 5194 

An absorbed power-law model is a poor fit to the PN 
and MOS spectra of the nucleus of NGC 5194 (C = 
4117.2 for 409 degrees of freedom (dof), GOF=100.0%). 
To show the significant spectral features, we fitted the ab- 
sorbed power-law model in the 3 — 5 keV and 7—10 keV 
bands and plotted the ratio of the observed data in the 
0.3 — 10 keV band and the best-fit power-law in FigureEl 
A strong soft X-ray excess below 1.5 keV and a narrow 
iron Ka line at 6.4 keV are evident in the spectrum. The 
X-ray image, shown in Fig. ^ clearly shows presence of 
extended emission which we interpret as originating from 
optically thin thermal plasma. Accordingly, we added a 
mekal component to describe the soft component and 
a Gaussian line model at 6.4 keV. The fit improves 
significantly (C = 452.3 for 404 dof, GOF = 88.5%), 
however, the soft X-ray emission is not well described. 
There are still line-like residuals below 2 keV, suggesting 
either the hot gas responsible for the soft excess emis- 
sion has a temperature structure and/or its abundance 
ratio is non-solar. Including an additional mekal compo- 
nent (model 1) further improves the fit {C = 414.7 for 
401 dof; 47.8%), and removes the soft residuals seen ear- 
lier. The best-fit temperatures are kT = 17lt^^ eV and 
kT — 610li[4 eV, and the best-fit abundances, relative 

to solar, are 0.20to:i4 and 0.18to:o2 the two mekal 
components. The observed data, the best-fit model and 
their ratio are plotted in Figure [T] We also tested the 
possibility of the non-solar abundance ratio by replacing 
the two mekal component by a vmekal component (model 
2) which allows the abundance ratios of individual ele- 
ments to vary. This model also improved the fit over that 
involving a single mekal component. The best- fit param- 
eters for both models are listed in Table|3| Model 1 and 2 
describe the data statistically equally well {GOF — 48%, 
and 47% for model 1 and 2, respectively). Both models 
result in an extremely flat hard X-ray slope (< 0.8), and 
strong iron Ka line (EW ^ 3 keV) from neutral material. 

The flat spectrum and strong iron line strongly sug- 
gest that the hard X-ray emission from the LLAGN 
in NGC 5194 is dominated by reflection from a cold 
material. Previous BeppoSAX observations of M 51 
showed the LLAGN is obviously Thomson thick {Nh ~ 
5.6tt °6 X 10^'' cm-2; Fukazawa et al. 2001). There- 
fore, we replaced the simple power-law model with the 
Compton reflection model pexrav (Magdziarz & Zdziarski 
1995). This model calculates the reflected spectrum from 
a neutral disk exposed to an exponentially cutoff power- 
law spectrum. Since XMM-Newton does not cover the 
expected peak (20-40 keV) of the reflection component. 
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it is not possible to tightly constrain all the parameters 
of the reflection model using this data alone. Instead, 
we fix the power-law photon index at 1.9, cutoff energy 
of the primary power law at 200 keV, disk inclination at 
70 degrees and the abundance of heavy elements at 0.1 
relative to solar value. We also fixed the relative amount 
of refiection, R= 1, corresponding to an isotropic source 
above a refiecting plane. The only variable parameter of 
the pexriv model was the normalization of the incident 
power law. We chose a high inclination angle due to the 
Thomson-thick nature of the LLAGN. We also set the re- 
flection model to produce the reflection component alone. 
This corresponds to a situation in which the primary X- 
ray emission is completely absorbed by a column density 
in excess of ~ 10^"* cm~^ along the line of sight and we 
observe the X-ray emission scattered by the matter near 
the nucleus. We also include a narrow Gaussian line 
model to describe the Fe K emission. Fitting the above 
model (model 3) resulted in no strong apparent residuals. 
Model 3 consisting of two mekal components, Compton 
reflection, and a narrow Gaussian resulted in C = 420.9 
for 400 dof and GOF = 61.8%. We consider this model 
to be the best-fit model as it is more physical than th ear- 
lier model involving simple power-law. The best-fit mekal 
temperatures are 177^11 eV, 614tn eV. The iron line 
equivalent width is 2.5^}'4 keV, similar to that measured 
from Chandra observations (Levenson et al. 2002). The 
best-fit parameters are listed in Table |21(Model 3). The 
observed flux is 8.2 x 10"^^ erg s~^ in the 0.3 - 10 keV 
band, and 2.0 x 10~^^ erg s~^ in the 2 — 10 keV band, 
which correspond to luminosities 6.9 x 10^^ erg s~^ in the 
0.3 — 10 keV band, and 1.7 x lO'^^ erg s~^, respectively. 

5.2. ULXs in M 51 

We performed joint fits for the PN and MOS spec- 
tra of each ULX with the following spectral models: 
(i) a power law (PL), (m) multicolor accretion disk 
blackbody (MCD), [Hi) thermal plasma (mekal), (iv) 
PL-fMCD, (v) PL+mekal. The models also incorporate 
photoelectric absorption [wabs with a minimum value of 
1.5 X 10^° cm~^, the Galactic absorption in the direction 
of M 51; Stark et al. 1992). We assume an abundance of 
0.1 relative to the solar in all the mekal components. The 
point spread function for the two sources 82 and 37 par- 
tially coincide with PN CCD chip gaps and/or bad pixel 
columns, prompting us to restrict ourselves to using only 
MOS data for the spectral analysis of these sources. One 
or more of the above models provided a reasonably good 
fit to all the ULXs except for sources 9, 26 and 69. The 
best-fit models and the unfolded X-ray spectra of these 
ULXs are shown in Figure |H1 The best-fit parameters 
are listed in Table 4. We have also listed the 0.5 — 8 keV 
observed flux and luminosity for the PN data in Table 4. 
The errors on flux were estimated at 90% level based on 
1000 sets of parameter values drawn from the distribu- 
tion that is assumed to be multivariate Gaussian centered 
on the best-fit parameters with sigmas from the covari- 
ance matrix. The fluxes measured with PN and MOS 
cameras are within the 90% errors quoted in Table 4. 
The mekal model alone does not provide a good flt to 
the spectra of all ULXs except source 12, therefore the 
best-fit parameters are not listed in Table 4. 

X-ray spectra of sources 5, 37 and 41 in NGC 5194 are 



equally well described by a PL-I-MCD or a Fh+mekal 
model, while source 63 prefers PL+mekal model. Two 
ULXs source 82 in NGC 5194 and source 12 in NGC 5195 
only requires a single spectral component. Source 82 is 
described by a simple power law and does not statistically 
require any soft component at a level > 99%. Source 12 
is well described either by a simple power law or a mekal 
model. In Table 4, we have listed the reduction in the 
C-statistic value (AC) for Ap additional parameters for 
the addition of a soft component {mekal, MCD or both, 
see below) to the simple power-law model. We have also 
listed the statistical significance derived from the maxi- 
mum likelihood ratio (MLR) test. The 99% criterian for 
selecting the more complex model is significance < 0.01. 
As a caveat we note that the MLR test is statistically not 
robust in cases that involve testing a hypothesis that is 
on the boundary of the parameter space e.g., the null val- 
ues of the normalizations of MCD, mekal and Gaussian 
lines occur in the simple power law model (Protassov et 
al. 2002). The results of the MLR test presented here 
should be treated only as indicative. 

The spectra of sources 9, 26 and 69 are not satisfacto- 
rily described by any of the four models mentioned above. 
Source 9 is the softest ULX in our sample. Although not 
a good fit, the PL-I-MCD model is marginally better than 
the PL-f mefcaZ model for this ULX (see Table 4). Adding 
a mekal component to the PL-f-MCD model improves the 
fit (AC = —9.7 for two additional parameters), and re- 
suhs in a good fit (C = 59.8 for 50 dof, GOF = 46.9%). 
Replacing the MCD component by another mekal com- 
ponent in the PL-fMCD-f mefca/ model slightly worsened 
the fit (C = 67.9 for 51 dof; GOF = 78.2%). 

The ULX source 26 is located in the region of strong, 
extended soft X-ray emission. Its X-ray spectrum is not 
well described either by a simple power law or an opti- 
cal thin plasma. Similar to the spectral modeling of the 
LLAGN, we used multiple mekal components to describe 
the soft emission. Two mekal components are statisti- 
cally required to explain the soft X-ray excess emission. 
The addition of the second mekal component is signifi- 
cant at a level > 99.99% (AC -28 for two additional 
parameters). In Figure |51 we show the ratio of the EPIC 
data and the best-fit model consisting of two mekal com- 
ponents and an absorbed power-law. An iron Ka line at 

6.4 keV is clearly seen in the spectrum. The addition 
of a Gaussian line improves the fit (AC = 19.2 for three 
additional parameters). The best-fit resulted in C = 278 
for 295 dof {GOF = 9.3%). The observed EPIC data, the 
best-fit model and their ratio are plotted in Figure ^| 
The best-fit parameters are listed in Table 

The above fit clearly shows the presence of a Gaus- 
sian line. The best-fit line parameters are Eii„e = 
e.SStoll keV, a = 158^^^^ cV, linefiux/i„,e = 2.3tl lx 
10^^ photons cm^^ s^^ and EW — 550 eV, consistent 
with an iron Ka line from neutral material. The LLAGN 
in NGC 5194, a Compton-thick AGN, also shows strong 
{EW ~ 3 keV; see also Levenson et al. 2002). The 
ULX NGC 5194#26 is separated by only 28" from the 
LLAGN. The fractional encircled energy for the PN is 
- 80% at a redius of 28". The location of NGC 5194#26 
within the wings of the point spread function and the 
very strong narrow iron line of LLAGN raise doubt if 
the iron Ka emission inferred from the ULX is sim- 



6 



ply cross-contamination from the LLAGN. To investigate 
this possibihty, we extracted a background spectrum us- 
ing circular regions with centers at the same distance 
from the LLAGN as that of NGC 5194#26. The radii 
of the two circles were chosen to be the same as that 
of the circular region used to extract the spectrum of 
NGC 5194#26. We used the PN data only and car- 
ried out the spectral analysis in the 2 — 10 keV band. 
A simple absorbed power-law model provided C = 83.8 
for 77 dof and GOF = 57.0%. Addition of a Gaus- 
sian line improved the fit (AC — —13.4 for three ad- 
ditional parameter). The best-fit parameters are F = 
2.4t^j, npL = 3.0l°j X 10-4 photons cm'^ s"! keV-^, 
Eune= 6.3t°;^ keV, a = 210^1% eV, fune = 2.7^11 x 
10"^ photons cm^^ s^^, EW — 690 eV. Replacing the 
Gaussian line with an accretion disk line (the diskline 
model in XSPEC) resulted in C = 69.5 for 74 dof. The 
best-fit parameters are -EHne — 6.31;q4 keV, disk in- 
clination angle, i = 26^26 degree, fune = 'i-Q-i^g x 
IQ-^ photons cm^^ s^^. The emissivity index, inner and 
outer radii were kept fixed at —2.5, 6rg and lOOOr^, 
respectively, where rg is the gravitational radius. The 
addition of the diskline to the simple power law model 
improved the fit at a significance level of 99.7%. We 
also tested if the broad iron line-like feature could be 
described in terms of an absorption edge. Addition of 
an edge at ^ 7.1 keV to the absorbed power-law model 
improved the fit (C = 70.6 for 75 dof). The best-fit pa- 
rameters are Nh = 2.5t|^ x 10^2 cm^^ F = 1.2j:°j, 

edge energy E^dge = 7.lt°ol keV and r = lA^-^^. Thus 
the 2 — 10 keV EPIC PN spectrum of source 26 is equally 
well described by an absorbed power law and a broad iron 
line at ~ 6.3 keV or an iron K-edge at ~ 7.1 keV. It is 
still possible that the background correction, we adopted 
here, may not be accurate due to poor signal-to-noise of 
the data, the presence of extended X-ray emission and 
moderate spatial resolution of XMM-Newton. Future 
long X-ray observations with Chandra and XMM-Newton 
will help detecting the line and/or edge unambiguously. 

Source 69 is located in a region of extended soft X- 
ray emission in an spiral arm. It is likely that X- 
ray emission from source 69 is contaminated by the 
extended emission. Adding another mekal component 
to the Ph+mekal model improved the fit significantly 
(AC = 52.7 for two additional parameters; C — 191.9 
for 200 dof; GOF = 18%). The best-fit parameters are 
listed in Table 4. Replacing one of the mekal compo- 
nent by the MCD component slightly worsened the fit 
(C = 200.9 for 199 dof, GOF = 29.2%). However, both 
the fits are acceptable. The spectral data and the best-fit 
model Ph+mekal+ mekal are plotted in Figure |S1 

We compared the soft excess flux in the 0.3 — 2 keV 
band, modeled as mekal component above a power-law, 
of ULXs with that of the diffuse emission in the sur- 
rounding regions. We estimated the flux of the diffuse 
emission from the count rates measured in annular re- 
gions centered at the source positions or nearby circular 
regions. We used the best-fit mekal temperature, esti- 
mated from the Plj+mekal model for an ULX, to con- 
vert the appropriately scaled count rate of the nearby 
diffuse emission into the flux in the 0.3 — 2 keV. For 
all ULXs except source 69, the nearby diffuse emission 



is comparable to the soft excess emission, suggesting 
that the soft excess emission may not be associated with 
these ULXs . The soft excess emission of source 69 
{f mekal ^ 6.6 X 10^"^ erg cm-2 s"^) is about a factor 
of four stronger than the surrounding diffuse emission in 
the 0.3 — 2 keV band. Thus, at least a part of the soft 
excess emission is intrinsic to the ULX. 

5.2.1. The comptonized accretion disk model 

The power law and the MCD model used to fit the ULX 
spectra above are only a mathematical approximation of 
the real physical spectra and do not provide much phys- 
ical insight. Therefore, we also tried a spectral model 
based on 'real' physics of an accretion disk-corona sys- 
tem. This model called as the 'comptonized accretion 
disk (CMCD)' has been used to describe the spectra of 
six bright ULXs (Wang et al. 2004). The CMCD model 
was constructed by Yao et al. (2003). The model as- 
sumes a thermal energy distribution for electrons in a 
spherical corona around an accretion disk. The parame- 
ters of the CMCD model are the radius of the spherical 
corona (Re), electron temperature (Tc), coronal optical 
depth (tc), the temperature of the inner accretion disk 
{kTin), and disk inclination angle (i). The model is im- 
plemented as a standard XSPEC table model (see Yao 
et al. 2003 and Wang et al. 2004 for more details). 

We fitted the CMCD model to the spectra of all the 
ULXs. It was not possible to constrain the parameters 
for all but one ULX source 26, due to poor signal-to-noise 
and/or best-fit parameters outside the tabulated ranges. 
Due to the contamination of the extended soft X-ray 
emission, source 26 statistically required two mekal com- 
ponents in addition to an absorbed CMCD model. The 
addition of the first and second mekal components im- 
proved the fit by AC = —491.1 and —19.1, respectively, 
for three additional parameters in each case. The mekal 
temperatures are similar to that derived earlier using the 
combination of two mekal and a power law model. The 
best-fit parameters are listed in Table For the best-fit 
CMCD model, the inner disk temperature is 291^g2 eV. 
The disk inclination angle is poorly constrained. 

6. INDIVIDUAL ULX CHARACTERISTICS 

In this section we describe the X-ray properties of in- 
dividual ULXs. 

6.1. NGC 5194 source 5 

This source is located in the western outer spi- 
ral arm of NGC 5194, coincident with the ROSAT 
HRI source 4 (Ehle et al. 1995). Chandra resolved 
this source into two discrete sources NGC 5194#5 
(CX0M51 J132939. 5-1-471244) and NGC 5194#6 
(CX0M51 J132940.0-H471237). NGC 5194#5 is 
brighter than NGC 5194#6 by factors of 1.5 and 11 
in the 0.5 — 8 keV and 2 — 8 keV bands, respectively 
(Terashima & Wilson 2004). Our XMM-Newton 
spectrum of this source requires a soft X-ray excess 
component that is equally well described either by an 
MCD or a mekal component. This soft excess component 
contributes ^ 42% to the total flux in the 0.5 — 2 keV. 
Therefore, the detection of soft X-ray excess emission 
above the power-law component in the XMM-Newton 
spectrum could be due to the contamination of the 
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softer NGC 5194#6 and the ULX spectrum is likely to 
be a simple power-law. 

6.2. NGC 5194 source 9 

This is the softest ULX in M 51 (see Fig. [TJ. This 
source is coincident to the ROSAT source 5 (Ehle 
et al. 1995) and Chandra source NGC 5194#9 which was 
a super-soft source (Terashima & Wilson 2004). Chan- 
dra detected no photons above 1 keV from NGC 5194 
source 9 and its spectrum was well described by an 
MCD model with kT ~ 100 eV. XMM-Newton EPIC 
MOS detected this source at 6a and 3cr levels above 
1 keV and 2 keV, respectively. The MOS spectra of this 
source statistically require three components: a power- 
law (r ~ 1.4), an MCD {kT - 130 eV), and a mekal 
{kT ~ 340 eV). Afternatively, the MOS spectra can 
also be described by a power law (F ~ 1.6) and two 
mekal components {kT ^ 80, 306 eV). The X-ray flux 
of source 5 was a factor of ~ 3 higher during our XMM- 
Newton observation than that during the Chandra obser- 
vations. The soft X-ray excess (MCD -I- mekal) flux is 
7.4 X 10^^^ erg cm^^ s^^ in the 0.5 — 8 keV band which 
is similar to the total flux observed by Chandra in June 
2001. 

6.3. NGC 5194 source 26 

This is the nearest ULX to the LLAGN in NGC 5194. 
It is located about 28" west of the LLAGN and its soft 
X-ray emission is heavily contaminated by extended soft 
emission. This source was not detected in the ROSAT 
HRI possibly due to the extended soft X-ray emission. 
This ULX is brighter than the LLAGN in the 1.5-10 keV 
band (see Fig. EJ. Chandra also detected this source as 
the hardest ULX in M 51 (Terashima & Wilson 2004). 
In the two Chandra observations, this source showed in- 
teresting spectral behavior. The first Chandra observa- 
tion revealed an absorbed power-law spectrum and emis- 
sion lines at ~ 1.8 keV, 3.24 keV, - 4.03 keV, and - 
6.65 keV, most likely the Ka lines of Si, Ar, Ca, and Fe, 
respectively (Terashima & Wilson 2004). In the second 
Chandra observation, no emission lines were clearly de- 
tected although the source was brighter by ~ 50%. The 
XMM-Newton spectra of this source are well described by 
an absorbed power-law (F ~ 2.4, Nh ~ 7.1 x 10^^ cm^^) 
and two mekal components {kT ^ 250 eV, ^ 590 eV). 
The absorbed power-law is consistent with that observed 
by the second Chandra observation. The two mekal com- 
ponents represent the contribution of the extended soft 
emission. We do not detect any emission lines except the 
iron line or an edge. Future deep X-ray observations will 
establish if the line is real. The EPIC FN spectrum of 
source 26, after appropriately correcting for the contribu- 
tion of extended emission, requires either a moderately 
broad {a ~ 200 eV) iron Ka line or a strong iron K edge 
at - 7.1 keV. 

6.4. NGC 5194 source 37 

This source is coincident with a southern spiral arm 
of NGC 5194 that displays extended UV emission (sec 
Fig-El- This source was not detected in the first Chandra 
observation but was detected in the second observation 
(Terashima & Wilson 2004). During our XMM-Newton 
observations, the ULX was ~ 50% brighter than dur- 



ing the second Chandra observation. The Chandra spec- 
trum was described by a simple power law, while the 
XMM-Newton spectra statistically require a soft compo- 
nent that is equally well described by a mekal or an MCD 
component, in addition to a power law. The power law is 
flatter (F = 0.8lg g) than that derived from the Chandra 

data (F = 1.55^q ^[5). It is likely that this discrepancy is 
attributable to the different number of spectral compo- 
nents used to describe the Chandra and XMM-Newton 
spectra. 

6.5. NGC 5194 source 4I 

This ULX is closer to the nucleus of NGC 5195 than 
to the nucleus of NGC 5194 but may be associated 
with the outer edge of an outermost spiral arm to the 
north of the nucleus of NGC 5194. This source was 
also detected in the ROSAT HRI (source 10 in Ehle 
et al. 1995). Chandra resolved this source into two 
sources NGC 5194#41 (CX0M51 J132953.7-t-471436) 
and NGC 5194#42 (CX0M51 J132953. 8+471432). 
NGC 5194#41 is a factor of '-^ 2 more luminous than 
NGC 5194#42. Remarkably, the XMM-Newton sonice is 
fainter by 50% than the Chandra source NGC 5194#41 
alone. Chandra spectra of this ULX are consistent with 
either a power-law or an MCD and no variability be- 
tween the two observations. XMM-Newton spectra sta- 
tistically require either an MCD {kT ^ 270 eV) or a 
mekal {kT ^ 340 eV) in addition to a power-law. Since 
the Chandra source NGC 5194^42 is an absorbed source, 
the thermal component inferred from the XMM-Newton 
data is likely a genuine component. 

6.6. NGC 5194 source 63 

This is faintest ULX in M 51 with a 0.5 - 8 keV lu- 
minosity of ^ 10^^ erg s~^. This source was fainter by a 
factor of ~ 2 in the 0.5 — 8 keV band during the Chandra 
observations and was not considered to be a ULX. We 
consider this source to be a ULX at its increased flux. 
It is possible that the source is a stellar mass black hole 
accreting at or near the Eddington rate. The ULX is 
located to the south-east region coincident with a spiral 
arm, displaying extended UV emission. The X-ray spec- 
trum of this source shows soft excess emission above a 
power law (F ~ 2.0) and is statistically better described 
by a mekal plasma than an MCD. The presence of an 
extended UV source and optically thin thermal plasma 
both suggest that this ULX may be associated with a re- 
gion of strong star-formation or unusually X-ray bright 
supernova remnant. 

6.7. NGC 5194 source 69 

This ULX is well known for the tentative detection of 
a 2.1 hr periodicity in its X-ray flux (Liu et al. 2002). 
It is located in a region of extended soft X-ray emission 
and strong UV emission in an spiral arm to the north- 
east of the nucleus. The ULX is closer to the nucleus of 
NGC 5195 than to the nucleus of NGC 5194 but is clearly 
associated with NGC 5194. This source was detected in 
the ROSAT ma (Ehle et al. 1995). The ULX showed 
drastic long term X-ray variability, its 0.5 — 8 keV lumi- 
nosity declining by a factor of ~ 75 from the flrst to the 
second Chandra observation. During our XMM-Newton 
observation, the ULX had brightened by a factor of ^ 50 
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in the 0.5 — 8 keV band. Source 69 also shows variabihty 
on short-times scales. The first Chandra observation ten- 
tatively detected two peaks separated by ~ 7000 s. Liu 
et al. (2002) interpreted these features as a 2.1 hr peri- 
odicity. The number of counts was too low to detect any 
such periodicity during the second Chandra observation. 
During a bright flux state during the XMM-Newton ob- 
servation, the source again shows three possible peaks 
separated this time by ~ 6000 s. We detect a peri- 
odic signal at a significance level of 2a with a period of 
5925 ±200 s from the EPIC light curves (see section lTn 
for further discussion). The XMM-Newton spectrum of 
this source is well described by either of the following 
three component models: (i) a power law (F 1.2) and 
two mekal components {kT ~ 180 eV, 680 eV); (ii) a 
power law (F ~ 1.2), an MCD (kT ~ 170 eV), and a 
mekal (kT ~ 690 eV). 

6.8. NGC 5194 source 82 

This is the brightest ULX in M 51 with an observed 
luminosity of - 2.2 x 10^^ erg s'^ in the 0.5 - 8 keV 
band. The ULX is coincident with an outer spiral arm 
east of the nucleus. This source was previously detected 
with Einstein HRI (Palumbo et al. 1985), ROSATRRl 
(source 14 in Ehle et al. 1995), BeppoSAX (Fukazawa et 
al. 2001), and Chandra (Terashima & Wilson 2004). The 
source had dimmed in the second Chandra observation 
but was brighter by 35% in the XMM-Newton observa- 
tion, a level similar to that seen during the first Chandra 
observation. Both Chandra and XMM-Newton spectra 
of this ULX are consistent with an absorbed power law. 

6.9. NGC 5195 source 12 

This is the only ULX in NGC 5195 that is weU sepa- 
rated from the nucleus. The ULX is ~ 70" away to east 
of the nucleus. The ULX was earlier detected with the 
ROSAT HRI (source 13 in Ehle et al. 1995). Source 
12 is the only ULX in M 51 to have X-ray spectrum 
that is well described by a single spectral component ei- 
ther a simple power-law (F ^ 1.6 ± 0.2) or mekal plasma 
{kT r^7.7 keV). The ULX is not resolved spatially with 
Chandra. Modeling the spectrum as an optically thin 
plasma requires an unrealistically high temperature and, 
as a result, we prefer the non-thermal power-law inter- 
pretation despite the equally good spectral fits. The no- 
table lack of an optical/UV counterpart for the source 
strongly suggests that this source is not a background 
AGN, particularly given its low intrinsic X-ray absorp- 
tion {Nh ~ 5.6 X 10^° cm-2) 

7. DISCUSSION 

We have performed the temporal and spectral analysis 
of the LLAGN and nine ULXs in M 51 using a XMM- 
Newton observation. Below we discuss the main results 
obtained. 

7.1. The LLAGN of NGC 5194 

The X-ray spectrum of the LLAGN is a typical of 
Thomson-thick Seyfert 2 galaxies, namely an extremely 
flat continuum, a strong and narrow line from neutral 
iron, and unabsorbed soft X-ray emission. Our spectral 
modeling shows that the fiat continuum can be described 
as the reflection of the primary power-law (F = 1.9) by 



cold and dense material, assuming the cold material sub- 
tends a 2tt solid angle to the central source, the lumi- 
nosity of the primary continuum is then estimated to 
be - 10"*° erg s'^ in the 0.3 - 10 keV band. This lu- 
minosity is about two orders of magnitude lower than 
that for Seyfert 1 galaxies, thus confirming the nucleus 
of NGC 5194 as an LLAGN as already established with 
previous BeppoSAX and optical observations. Since the 
LLAGN is nearly Thomson thick, the soft X-ray emission 
in the spectrum must be associated with the extended 
soft emission seen in Fig. ^ 

7.2. ULXs in M 51 

There are nine ultra- luminous X-ray sources in M 51 
with the observed luminosity of ^ 10^^ erg s~^ in the 
0.5 - 8 keV band. Three ULXs, source 9, 41 and 63, 
with Lx ~ 10'^^ erg s~-^ (0.5 — 8 keV) may actually be 
stellar mass black holes (~ 10 — 2OM0). The brightest 
ULX in M 51, source 82, has a 0.5 — 8 keV luminosity of 
2.2 X 10'^^ erg s~^, lower than the Eddington luminosity 
(~ 2.6 X 10^9 erg s'^) for a 20 M0 black hole. Thus, aU 
ULXs in M 51 could plausibly be stellar mass black holes 
at around (or slightly above) the Eddington limit. 

The larger collecting area of XMM-Newton compared 
to Chandra has enabled us to investigate the nature of 
ULXs in greater detail. While Chandra spectrum of in- 
dividual ULXs in M 51 is well described by a single spec- 
tral component, XMM-Newton spectrum of most ULXs 
in NGC 5194 require at least two spectral components. 

7.2.1. Cool accretion disks 

X-ray spectra of two ULXs, sources 9 and 69, require 
soft components, modeled by either mekal plus MCD or 
two mekal components. The MCD components suggest 
inner disk temperatures of kTin ~ 130 eV (source 9) and 
~ 170 eV (source 69). These temperatures are lower 
than that found for ULXs studied with ASCA (Mak- 
ishima et al. 2000) and Galactic X-ray binaries, but 
are similar to that found for ULXs studied with XMM- 
Newton (Miller et al. 2003; Miller, Fabian, & Miller 2004; 
Dewangan et al. 2004). The XMM-Newton observation 
of M 51 supports the increasing trend that ULXs have 
cooler accretion disks than GBH binaries. The tempera- 
ture of the innermost regions of an accretion disk is given 

by 

(1) 

(see e.g., Makishima et al. 2000). Here n 2± 1.7 is a 
spectral hardening factor. The factor f ~ 0.4 takes into 
account the fact that the radius of maximum tempera- 
ture is larger than the innermost stable orbit, and a — 1 
for a Schwarzschild geometry or 1/6 for a Kerr geometry. 
Thus objects with lower MCD temperatures must have 
either higher black hole masses or lower accretion rates. 

Assuming a relative accretion rate of ^0.1, equa- 

tion n implies black hole masses of ^ 73OOM0 (source 
9) and - 25OOM0 (source 69). 

The results of the CMCD model fit to the spectrum 
of source 26 suggests that this ULX has cool accretion 
disk {kT — 29lig2 ^V). Using the normalization and the 
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inclination angle of CMCD model, it is possible to derive 
a BH mass as detailed in Wang et al. (2004). Using this 
procedure, we derive a black hole mass of ~ 28OOM0 for 
the ULX sources 26. 

7.2.2. ULXs with a power-law spectrum 

The X-ray spectrum of two ULXs, source 82 and source 
12, can be described by a simple power-law models with 
a photon indices of ~ 2.4 and ^ 1.7, respectively. Source 
12 does not show any signature of emission from an ac- 
cretion disk. The power-law spectral form of this ULX is 
similar to that of XRBs in their low state (Terashoma & 
Wilson 2004) or type 1 AGNs (excluding the narrow-line 
Seyfert 1 galaxies). If the above analogy is correct, then 
the ULXs with simple power-law spectrum may imply 
very high mass (> a few times 1OOM0) black hole. A 
power-law spectrum may also result from beamed X-ray 
emission. The X-ray spectrum of source 82 is very steep 
and may be similar to that of Galactic XRBs in their very 
high state (see Terashima & Wilson 2004 and Roberts et 
al. 2004 for a detailed discussion on the power-law spec- 
tral form of ULXs) . It is also possible that ULXs with a 
simple power-law spectrum are background AGNs. 

7.2.3. Spectral variability 

It is not possible to compare directly the best-fit spec- 
tral form derived from the Chandra and XMM-Newton 
observations due to the different number of spectral com- 
ponents used to fit the spectra. Therefore, we compare 
the best-fit model derived from the Chandra, usually 
a single component model, and the corresponding sin- 
gle component model best-fit to the XMM-Newton data. 
The single component model is not necessarily the best- 
fit model to the XMM-Newton data. FigurelTTl compares 
the spectral parameters illustrating the spectral variabil- 
ity of ULXs in M 51. The steeper spectrum of source 5 
during the XMM-Newton observation may be due to the 
contribution of the soft source NGC 5194#6 detected 
with Chandra (see section 6.1). Source 9 shows a very 
soft X-ray spectrum and is consistent with no spectral 
shape variability. The detection of the power law com- 
ponent in addition to the soft component in the spectra 
of source 9 may be due to the larger collecting area of 
XMM-Newton particularly at higher energies. Sources 
26, 37, 63 and 12 do not show significant variability in 
their spectral shape. Sources 41 and 82 show slightly 
steeper X-ray spectra compared with the earlier Chan- 
dra observation. 

Source 69 continues to show drastic spectral and fiux 
variability. The source was in a high/hard state (power 
law r = 1.24tE;j^ jy^cD = 2.3^0^) in June 
2000 and subsequently made a transition to an ex- 
tremely soft /low state (power law F > 5.1 or MCD 
kT = 0.17;^;J^) in June 2001 (Terashima & Wilson 
2004) . The ULX made another transition back to a state 
similar to its earlier high/hard state, in January 2003. 
The spectral transitions of source 69 are completely dif- 
ferent from the high/soft to low/hard spectral transitions 
usually observed from Galactic XRBs but similar to those 
observed from some ULXs (The Antennae, Fabbiano et 
al. 2003; Holmberg II X-1, Dewangan et al. 2004; an 
ULX in NGC 7714, Soria & Motch 2004). The ULXs 
that show distinct spectral transitions may a comprise a 



distinct class of ULXs. The physics of these transitions 
is yet to be understood clearly. 

7.2.4. Periodicity 

Source 69 is the only ULX in M 51 that shows indica- 
tions for a possible periodic variations in its X-ray flux. 
The periods detected from our XMM-Newton observation 
(P = 5925 ± 200 s at a 2cr level) and the first Chandra 
observation (P — 7620 ±500 s) are significantly different. 
The two observations are separated by ~ 1.5 yrs. The 
variation in period, if real, clearly argues against orbital 
periodicity, and suggests the presence of quasi-periodic 
oscillations. Liu et al. (2002) discussed the possibility 
that the observed periodicity is due to orbital variations 
and as noted by Liu et al. (2002), only low mass X-ray 
binaries (LMXBs) can have orbital periods as small as 
2 hr with high luminosities. In this case Liu et al. (2002) 
estimate the mass of the donor to be ~ O.23M0. The 
spectrum of source 69 derived from the first Chandra ob- 
servation was found to be consistent with an absorbed 
multicolor disk model with inner disk temperature of 
kT = 2.101q4q -'^^V ^^'^ showed marginal improvement 
in the fit after including a single temperature blackbody 
component. Based on these results, Liu et al. (2002) do 
not rule out the possibility that source 69 is a neutron 
star. Our XMM-Newton spectra require three compo- 
nents - a power-law (F ^ 1.2), a mekal {kT ~ 680 eV) 
and another mekal {kT ~ 180 eV) or a cool accretion 
disk [kTin ^ 170 eV). The stellar mass black hole sys- 
tems and high luminosity neutron star binary systems, 
both have high disk temperatures {kT = 1 — 2 keV). The 
possible variation in the period, due to a QPO and the 
X-ray spectrum both are inconsistent with the LMXB 
picture. Moreover, source 69 resides in a region that is 
likely to be undergoing active star formation, making it 
likely that the ULX is a high mass X-ray binary (HMXB) . 
The spectral transition of source 69 from a high-hard to 
low-soft state also is in disagreement with the behavior 
of most Galactic LMXBs (Terashima & Wilson 2004). 

If we interpret the variation of the apparent peri- 
odicity as resulting from QPOs, then QPO frequency 
of ~ Q.lbmHz is much lower than QTmHz QPO ob- 
served with RXTE from the Galactic black hole can- 
didate GRS 1915-1-105, which has a dynamically mea- 
sured mass of 14 ± 4M0 (Morgan, Remillard & Grener 
1997). The only ULX known to show QPO is M 82 X-1 
(Strohmayer & Mushotzky 2003). Its QPO frequency of 
MmHz is also lower than that of GRS 1915-f 105 but 
much larger than the implied frequency for source 69. It 
is not necessary the case that QPO frequency scales with 
the mass of the black hole, only twin peak QPOs with a 
frequency ratio of 3:2 have been established to scale with 
black hole mass (McClintock & Remillard 2003). The 
micro-quasar GRS 1915-1-105 is known to show QPOs at 
kHz as well as mHz frequencies, the lowest frequency of 
QPO observed from GRS 1915-H05 is 1.6 mHz (Morgan 
et al. 1997). Therefore it is not entirely unlikely the 
possibility of a ^ 0.15 mHz QPO from source 69. Long 
observations are required to clarify the presence of peri- 
odicity or QPO in this ULX. 

8. CONCLUSIONS 

We have analyzed XMM-Newton EPIC observation of 
M 51 performed in January 2003. Our main findings are 
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as follows: 

• We detected 9 ULXs, an LLAGN and extended soft 
thermal emission from M 51. One of the sources 
(source 63) is a new ULX, not identified earlier due 
to its previously faint flux levels. 

• All eight ULXs in NGC 5194 are located in or near 
spiral arms. Four ULXs are located in the regions 
of strong UV emission in the spiral arms, suggest- 
ing association of these ULXs with current star for- 
mation. 

• Two ULXs (sources 26 and 69) show evidence 
for short-term variability. Source 69 is the most 
rapidly variable ULX in M 51 and shows a pos- 
sible period of ~ 5925 s. All the ULXs except 
NGC 5194#12, show long-term variability by a fac- 
tors of few to several. 

• Most ULXs in M 51 show soft X-ray excess emis- 
sion, above a hard power law, that is modeled by 
a mekal plasma and/or an MOD component. How- 
ever the strength of the soft excess emission in- 
ferred from an ULX is comparable to the diffuse 
emission around individual ULXs excluding source 
69. Thus the soft excess emission is unlikely to be 
physically associated with any of the ULXs except 
source 69. The soft excess emission from source 69 
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is about a factor four stronger than the surrounding 
diffuse emission and is likely intrinsic to the ULX. 

• There is an indication of either a moderately broad 
(cr = 100 eV) iron Ka line or an iron K absorption 
edge both due to neutral iron in the spectrum of 
source 26. 

• The X-ray spectrum of LLAGN in NGC 5194 is 
extremely flat (F ~ 0.7) and shows a strong (EW~ 
2.5 keV), narrow, iron line from neutral material. 
The X-ray emission from the LLAGN in NGC 5194 
is dominated by reflection from cold material. 
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Fig. 1. — Left: Three color image of M 51 created from adoptively smooth XMM-Newton EPIC PN images. The three colors red, green, 
and blue correspond to the 0.2 — 0.7 keV, 0.7 — 2 keV, and 2 — 10 keV energy bands. North is towards the top and east is to the left of the 
image. The extended source at the top of the image is the companion galaxy NGC 5195. The bright extended source near the center of 
the image is the center of the galaxy NGC 5194 that consists of extended soft emission and two hard X-ray sources: the LLAGN and an 
ULX NGC 5194 source 26. Right: The EPIC PN image of M 51 in the 2-10 keV band showing the positions of 9 ULXs and the LLAGN. 




Fig. 2. — XMM-Newton ultra-violet image overlaid with the 1.5 — 10 keV band X-ray contours, derived from the PN image, showing the 
discrete hard X-ray sources. The contours are drawn at levels of 0.4, 0.7, 1.3, 1.9, and 2.6 counts pixel"^ (1 pixel = 0.953 arcsec). The 
image size 11.1' X 11.1'. North is up and east is left. Two central hard X-ray sources, about 20" apart, are clearly resolved by the EPIC 
instruments. The center of UV image coincides with the fainter of the two hard X-ray sources in the central regions. This fainter source is 
the LLAGN in NGC 5194. The brighter source is the ULX NGC 5194 source 26. 
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(a) Source 26 



(b) Source 69 



5000 



10* 

Relative time (s) 



1.5x10* 



Fig. 3. — Short-term EPIC PN light curves of ULXs in M 51 (a) source 26, and (b) source 69. The bin sizes are 520 s. The light curves 
have been corrected for the background contribution. 




Period (s) Phase 
Fig. 4.— Left: Periodogram of NGC 5194 source 69. Right: Folded XMM-Newton light curve of NGC 5194 source 69. 
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Fig. 5. — Long-term variability of ULXs in M 51 over a baseline of two and half years. The fluxes for the years 2000 and 2001 were 
derived from the two Chandra observations by Terashima & Wilson (2004). All but two ULXs 41 and 12 clearly varied. 
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Fig. 6. — Ratio of the observed EPIC data and the best-fit power-law model obtained by fitting the data in the 3 — 6 keV and 7 — 10 keV 
bands. 
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Energy (keV) 

Fig. 7. - The observed EPIC PN and MOS spectra of the nucleus of NGC 5194 and the best-fitting spectral model which a combination 
of two mekal plasma components, power-law, and a narrow Gaussian line at 6.4 keV and multiplied by the Galactic absorption. 
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Unfolded EPIC PN or MOS spectra of ULXs in M51. 
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Energy (keV) 

Fig. 9. — The ratio of the observed EPIC data and the best-fitting model consisting of three mekal components and a power law modified 
by the Galactic absorption. An iron Ka line at 6.4 keV is evident. PN, MOSl, and MOS2 data are marked by filled circles, open circles, 
and open squares, respectively. 
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Fic;. 10.— The observed EPIC PN and MOS spectra of the ULX NGC 5191#26 and the best-fitting spectral model - a combination of 
two mekal plasma components, a Gaussian line at ~ 6.4 keV, and an absorbed power-law. 
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Fig. 11. — Long-term spectral shape variability of ULXs in M 51. The best- fit power-law models (for all ULXs but NGC 5194 somce 9) 
and MOD model (for NGC 5194 source 9) to the Chandra data of June 2000 and June 2001, taken from Terashima &: Wilson (2004), and 
to the XMM-Newton data of January 2003 have been used. 
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sources fall on the PN chip gaps. The coordinate was derived from the full band image. 



TABLE 2 

Results of short-term variability tests. 
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TABLE 3 

Results of the joint spectral fittings to the EPIC PN and MOS spectra of M 51 nucleus 
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''(f) indicates a fixed parameter. 

''Model 1 is a combination of two mekal, a power-law, a narrow Gaussian multiplied the Galactic absorption model. Model 2 is the same 
as model 1 except that the two mekal components have been replaced by a vmekal component. Model 3 is the same as model 1 with the 
power-law replaced by a neutral reflection model pexrav. 

<=In units of 10"" / UeUHdV, where d is the angular size distance to the source ( cm), rie is the electron density ( cm ^), and uh 
is the hydrogen density ( cm"^). 

''Line flux in units of 10"® photons cm"^ s"'. 

"Observed flux in units of 10"'^ erg cm"^ s"' 
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TABLE 4 

Best-fit spectral model parameters derived from the joint PN/MOS data for the ULXs in M 51. The fluxes measured with 
the pn and mos data are similar within errors, and the quoted fluxes are from the pn data only. 
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40.970 


—53.9/4 


c c in — 11 

0.5 X lU 




r 


1.57(t) 






o'J_44,3Ub_!_39 


1 n+0.3 
1-0-0.2 


n a+0.2 
0-8_o.2 


67.9/51 


78.2% 


—45.8/3 


c: w in — 10 


37^ 


A 


1.57(f) 


1 -7 + 0.1 










104.9/61 


99.9% 








B 


1.57(f) 




1253_i3i 








z^2.o/d1 


1 nn nQ/ 

lOU.O/o 










in ^ + 15.9 


n a+0-4 
0-8_o.s 


220_67 




9 Q+0.4 


1 Q+0.3 
1-9-0.9 


bo.b/oo 


c:'7 ntV 
57.970 


—36.3/3 


a a \y in — 8 
D.5 X iU 




JJ 


1.57(t) 


n a+0.3 
0-8-0.5 




f.Qi+170 
691_94 


9 Q+0.5 

2-3_i.6 


1 Q+O.S 
1-9-1.2 


65.2/59 


65.2 


—39.7/2 


o A in — 9 
2.4 X lU 


41 


A 


1.57(f) 


1 o+O.l 

1-8-0.1 










169.1/129 


96.3% 








r> 


1.57(t) 




-1 1 -34 +M 

1134_g2 








Qno T / 1 on 

v>y2. / / izy 


1 nn nO/ 
100.0/0 








C 


A n+11.2 


1 n+O-3 
1-2-0.3 


^'^-100 




1 1+0.2 
1-1-0.5 


1 n+0.2 
I-OI0.4 


143.8/126 


59.8% 


—25.3/3 


1.3 X 10 




D 


1.57(f) 


1 r.+0-l 
-^-^-0.2 




0^1+86 
34II73 


1 1+0.1 
1-1-0.2 


n o+O.l 
0-9Io.2 


1 AC n 1 A OT 

146.6/127 


69.6% 


—22.5/2 


1.3 X 10~^ 


63 


A 




O q + O.l 

'^-•^-0.2 










137.0/101 


97.0% 








rS 


1.57(t) 




-700 + '^'^ 








QVfi c /I no 
( 0.0/ iUz 


1 nn no/ 
iUU.U /o 








C 


O n + 10.0 

3-9-3.9 


1 c + C..") 
-'--''-0.4 


no-j+lOO 
2»3_io5 




1 1+0.2 
1-1^5 


n n+0.2 
O-9I0.4 


116.3/99 


65.7% 


—20.3/2 


3.9 X 10~^ 




D 


l-57(/) 


1 Q+0.2 
-"--^-0.2 




6081 149 


1 1+0.2 
1-1-0.2 


n n+Ol 
0-9Io.2 


105.1/ 100 


42.01)70 


—31.9/1 


1.6 X 10 




A 
A 


1.67(t) 


9 1+0.1 










04Z.U/ 


1 nn nc% 
lUU.U/o 










1.57(1) 




OQ1 +42 
8yi_3g 








1 riQi Q /on/1 
iUyi.o/zU4l: 


1 nn nQ^ 
lUU.U/o 










1 n fi+7.0 


1 n + 0.2 
1-3-0.2 


IQOtgf 




9 Q+0.3 
^-"^-0.8 


1 n+0.2 
I-9I0.7 


oon n /oni 

220.0/201 


67.3% 


inn n /q 
— 122.0/0 


n o >^ in — 26 

2.8 X 10 






1 57('f'l 






368+52 




1 8+S-? 

^-"-0.2 


244.6/202 


87.6% 


— 106.4/2 


7.8 X 10~24 




E 


r- n+8.8 
°-^-5.2 


1 9+0.2 
-■--^-0.2 




689ti'^'' 


9 9+0.3 
^-^-0.6 


1 9+0-2 


200.9/199 


29.2% 


-141.1/5 


1.0 X 10^28 




F 


1.57(f) 


1 9+0.2 




18lli^, 6761112 


9 Q+0.3 


1 n+0.2 
^-^-0.4 


191.9/200 


17.1% 


-150.1/4 


1.9 X 10-31 




A 




9 4+0.2 






^-"-0.4 


9 9+0.4 


98.3/100 


32.4% 








B 


1.57(f) 




809+^^ 








164.2.3/101 


100.0% 








C 




1 4+0.6 


406l« 




9 7+0.1 
^- ' -2.0 




91.5/98 


12.8% 


-6.8/2 


0.033 




D 




9 9 + 0.1 
^-^-0.2 




'7Q'7+360 
-258 


^-"-0.6 


9 9 + 0.3 
^-^-0.5 


93.7/98 


21.3% 


-4.6/2 


0.100 


12 


A 




^-"-0.1 






-^-^-0.3 


1 o+"-2 


88.4/93 


21.3% 








B 


1.57(f) 




129011^° 








121.5/98 


92.2% 








G 


1.57(f) 




793012923 


-1 Q+0.2 


1 l+°-2 


87.7/94 


18.7% 







'»A= PL, B = MCD, C = PL+MCD, D = PL+mekal, E = PL+MCD+mefca/, F=PL+mekal + mekal, G=mekal. 
''Total absorption column density. The Galactic column is 1.5 X 10^0 cm" ^. 

'^Only MOS data were used for the spectral fitting. The source photons fall at the chip gap in the PN camera. 
"^Observed flux in units of 10—13 erg cm^"^ s^i and in the energy band of 0.5 — 8 keV band. 
"Observed luminosity in units of 1039 erg s"! and in the energy band of 0.5 — 8 keV band. 
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TABLE 5 

Results of the joint spectral fittings to the EPIC PN and MOS spectra of the ULX NGC 5194 source 26 



Component 


Parameter^ 


model 1'' 




model 2^ 


mekal 


kT{ eV) 






OKI 




abundance (solar) 


0.1(/) 




o-i(/) 




norm'^ 


i-8l J:^ X 10- 


-4 




mekal 


kT( cV) 


589140 




589139 




abundance (solar) 


o-i(/) 




0.1 




norm'^ 


3.41" 2 X 10 


-4 


3.41^ 2 X 10 


absorption 


TIT 1 1^99 9\ 

Nfi (xlO'''' cm ^) 






1 n P- + 1 4 
10-5I4.2 


power law 


r 


„ . ^4-0 39 
^■'±O_0.20 








norm 


3.9+f g X 10" 


-4 




Gaussian 


E kcV 






f, oo + O.ll 




cr ( eV j 






170_i36 




c 

J 




-6 


o 'vf" 1 ■ 5 -.J in — 6 

2.7 _^ 3 X lU " 




J?WT ( c,\T\ 

iIjVV ev^ 






DiO 




Klin [ eVj 






nni ~l"13 




kTc ( keV) 






cn n~l~35.T 

50-0Ii8.i 










6.65+|-2 




T 






Z.U_j 2 




id 










K" 






14*^ Q+383.9 


Total 


/x(xlO-" erg cm-2 s-l)s 






4 47+3.20 




LxixW^^ erg s-l)g 






3 76+2.69 
'°-1.96 




C/dof 


278.5/295 




275.0/291 




Goodness of fit 


9.3% 




10.6% 



''(f) indicates a fixed parameter. 

''Model 1 is a combination of two mekal, a power-law, a narrow Gaussian multiplied by an absorption model. Model 2 is same as model 
1 except that the power-law has been replaced with the comptonized MCD model. 
"In units of lO"" X 10"^*^;^ / nenndV, where d is the angular size distance to the source ( cm), Ue is the electron density ( cm 

and riH is the hydrogen density ( cm^^) 
"^Parameter pegged at both limits 

''CMCD normalization, K = (-^jU^i 
' ^ D/IO kpc ' 



^Observed flux and luminosity in the 0.5 — 8 keV band. 



